Quasiclassical description of a superconductor with a spin density wave 
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We derive equations for the quasiclassical Green's functions g within a simple model of a two-band 
superconductor with a spin-density- wave (SDW). The elements of the matrix g are the retarded, 
advanced, and Keldysh functions each of which is an 8 x 8 matrix in the Gor'kov-Nambu, the 
spin and the band space. In equilibrium, these equations are a generalization of the Eilenberger 
equation. On the basis of the derived equations we analyze the Knight shift, the proximity and 
the dc Josephson effects in the superconductors under consideration. The Knight shift is shown 
to depend on the orientation of the external magnetic field with respect to the direction of the 
vector of the magnetization of the SDW. The proximity effect is analyzed for an interface between a 
superconductor with the SDW and a normal metal. The function describing both superconducting 
and magnetic correlations is shown to penetrate the normal metal or a metal with the SDW due 
to the proximity effect. The dc Josephson current in an Ssdw / N / Ssdw junction is also calculated 
as a function of the phase difference ip. It is shown that in our model the Josephson current does 
not depend on the mutual orientation of the magnetic moments in the superconductors Ssdw and 
is proportional to sirup. The dissipationless spin current jsp depends on the angle a between the 
magnetization vectors in the same way [jsp ~ sin a) and is not zero above the superconducting 
transition temperature. 

PACS numbers: 74.20 Fg, 74.20 Rp, 74.72.-h, 75.40 Gb, 74.25 Dw, 74.25 Ha 



I. INTRODUCTION 

Coexistence of two or more order parameters in solids 
is one of the most intriguing phenomena in condensed 
matter physics. There are many systems where the or- 
dered phases are even antagonistic to each other. Nowa- 
days, a very popular subject of research are compounds 
with coexisting superconducting (SC) and magnetic (M) 
order parameters. 

It is already well known, that the exchange field re- 
sponsible for ferromagnetism destroys singlet Cooper 
pairsi but can generate^ triplet ones (odd in frequency 
and symmetric in space) that survives even in the pres- 
ence of a strong exchange field and strong impurity scat- 
tering. The theoretical prediction has been recently con- 
firmed experimentallyt£. 

In contrast, a magnetic order of the antiferromagnetic 
type can coexist with singlet Cooper pairs because the 
average magnetic moment (on distances of order of the 
size of the Cooper pairs) is close to zero. This case is real- 
ized in superconductors with a spiral magnetic structure 
(for a review, se&i^ and references therein). 

Another example of this coexistence are superconduc- 
tors with a spin-density wave (SDW)^^"^^. This type of 
superconductivity coexisting with the SDW is realized in 
quasi-one-dimensional conductors (for a review, seei^). A 
lot of attention is payed now to quasi-two-dimensional su- 
perconducting compounds - the so-called Fe-based pnic- 
tides ~ discovered recentlj*i^"— . Some of these Fe-based 
superconductors have rather high critical temperature Tc 
of the superconducting transition (up to 56 K observed 
in Gdi-xThxFeAsO at the doping level x = 0.3iS) and 
have much in common with the high cuprates (see, for 



example^Sr— and references therein). It has been estab- 
lished, both theoretically and experimentally, that at a 
certain critical temperature Tm (usually Tm > Tc) the 
compounds undergo a magnetic transition leading to the 
formation of a SDW. In a certain interval of the temper- 
ature and doping level x the SC and M order parameters 
may coexist2£iM"i21. 

The SDW in pnictides arises due to the nesting of the 
electron and hole pockets^'^>2^'^>2^. Each pocket has 
its own SC order parameter, Agj resp. A^, but the most 
energetically favorable state corresponds to the so-called 

S-i -pairing, characterized by the opposite signs of the 

SC order parameters: A^i = -An^^^^^^^L^liSS.. 

Microscopic theories for the electronic states in ma- 
terials with the superconducting and antiferromagnetic 
order parameters are mainly based on the equations for 
the Green's functions G which include both normal and 
anomalous ones corresponding to the SC and M order 
parameterSiS^"— These equations are written in the 
mean field approximation in an analogy with the or- 
dinary superconductors without specifying the micro- 
scopic mechanism of superconductivity and were solved 
under some approximations mostly in spatially homoge- 
nous case. On the other hand, several research groups 
have begun to study pnictides in nonhomogeneous struc- 
tures. For examples, the results concerning Josephson ef- 
fects in junctions based on pnictides have been published 
recently'^^"'^'^ . The Josephson effect in tunnel S/I/S junc- 
tions (/ denotes an insulator) consisting of the of multi- 
band superconductors can be investigated on the basis of 
the tunnel Hamiltonian and the equations for the Green's 
functions G^"— . However, this approach is not applica- 
ble to other types of the Josephson junctions. 
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In order to tackle the problems in nonhomogeneous 
cases it is much more convenient to employ the widely 
used quasiclassical approach based on using the so-called 
quasiclassical Green's functions g. These functions are 
obtained from the full Green's functions G by integration 
over the modulus of the momentum p in the vicinity of 
the Fermi surface. 

The quasiclassical Green's functions g are used in sit- 
uations when the parameters characterizing the system 
vary on the distances exceeding the Fermi wave length 
\p. The method of quasiclassical Green's functions has 
been developed in the theory of superconductivity^""— 
and turned out to be the most powerful and effective 
tool in dealing with the problems in nonuniform cases: 
vortices in superconductors, proximity effect in super- 
conductor/normal metal structures, etc. The equations 
for the quasiclassical Green's functions can be general- 
ized to the case of the charge-density-wave (CDW)^i^ 
and two-band superconductors^""—. They are very ef- 
ficient for describing superconductors in contact with 
ferromagnets^. 

At the same time, corresponding equations for the 
case of a two-band superconductor with the SDW are 
still lacking, although such equations can be derived 
under certain restrictions. Deriving these equations is 
very important because they can serve for the descrip- 
tion of the superconductors like Fe-based pnictides. The 
quasiclassical approximation is well justified for describ- 
ing the pnictides because the correlation lengths for 
both the superconducting and the magnetic correlations 
Cs,M ~ vh/ /S.,vh/WM are much longer than the Fermi 
wave length, where Wm Tm is a characteristic energy 
related to the magnetic order. 

In this paper, we derive the equations for the quasiclas- 
sical matrix Green's function g that describe a two-band 
superconductor with the SDW. These equations can be 
applied both to equilibrium and nonequilibrium states in 
homogeneous and nonhomogeneous cases and describe 
a broad class of phenomena in superconductors like Fe- 
based pnictides. We employ these equations to analyse 
the Knight shift, the proximity and the dc Josephson ef- 
fect in such superconductors. 

The paper is organized as follows. In Sec. II we write 
the Hamiltonian of the system in the mean-field approxi- 
mation in terms of the operators suitable for deriving the 
equations for the quasiclassical Green's functions. These 
equations are derived in Sec. Ill in the ballistic limit. 
We also present formulas for the observable quantities 
such as the SC and M order parameters (A and Wm), 
as well as for the charge and spin current density. In 
Sec. IV we study the Knight shift, i.e. the shift of the 
NMR (nuclear magnetic resonance) frequency due to the 



spin polarization of the s-electrons, for different orienta- 
tion of the external magnetic field with respect to the 
direction of the magnetization m in the SDW. In Sec. V 
the proximity effect will be analyzed in the vicinity of the 
interface between a superconductor with the SDW and 
a nonsuperconducting material (with or without SDW). 
Using a simple model, in Sec. VI we calculate the dc 
Josephson (jj) and the spin (j^p) current in the Joseph- 
son SsDw/N/SsDW junction, where Ssdw denotes a 
superconductor with an SDW, and A'^ is a normal metal. 
The obtained results will be discussed in Sec. VII. 



II. MODEL AND BASIC EQUATIONS 

We consider a simple model of a two band supercon- 
ductor with such a Fermi surface that not only the su- 
perconducting but also an SDW pairing is possible. The 
SDW pairing may originate from the nesting of certain 
parts of the Fermi surface and we assume that such parts 
exist. 

In such a situation one can have logarithmic contribu- 
tion not only from the Cooper channel but also from the 
particle-hole one. Solving this problem microscopically is 
not easy because one should perform rather complicated 
renormalization group calculations in order to get the in- 
formation about non-trivial phases at low temperature. 

We do not intend to discuss here microscopic mecha- 
nisms of the superconductivity and its competition with 
the SDW. Our goal is more modest: assuming that 
the superconductivity and the SDW coexist we consider 
them in the mean field approximation. In principle, one 
could obtain physical quantities of interest using the nor- 
mal and anomalous electron Green's functions written 
in the presence of the superconducting and SDW order 
parameters2^i21. 

Unfortunately, if the system is not homogeneous, solv- 
ing the equations for the electron Green's functions is not 
easy and therefore we develop the formalism of quasiclas- 
sical Green's functions g taking into account both the 
superconducting and the SDW order parameters. These 
equations are obtained as a result of a certain simplifi- 
cation of the original equations for the electron Green's 
functions. This approach is valid in the situations when 
physical quantities vary slowly on the distances of the 
order of the electron wavelength. 

To be close to the experimental results on pnictides we 
assume that the s-wave superconducting pairing inside 
the bands is most important. The spins of the electrons 
interact with both the exchange field of the SDW and 
the external field H. 

The Hamiltonian H of the system under consideration 
can be written in the form (see Refs i^^'^^ ) 
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W^zi [(0-3)0/3 cos6' + iai)at3 sine]c\^cip - T^Z2 [(0-3)0^ cos6' + iai)ap sm0]cl^C2p^ 



(2.1) 



where ^1,2 (p) is the kinetic energy in the 1 (hole) and 
2 (electron) bands, respectively, counted from the Fermi 
energy and A1.2 is the superconducting order parameter 
in these bands. 

The quantity Wmo = {Wm + W^m)/2 is the mag- 
netic order parameter describing the SDW (for simplicity 
we set the incommensurability wave vector q to zero), 
Wm = fJ-effn, where fJ-eff is the effective magnetic mo- 
ment of the free electrons participating in the formation 
of the SDW and m is the magnetization of the SDW. 

The terms Wzi = HiH and Wz2 = fJ-2H are the Zee- 
man energies in the presence of the external magnetic 
field with the components H = iJ(sin 0,0, cos 0), where 
fii, fi2 are the effective magnetic moments in the hole and 
electron bands, respectively. We introduce the Pauli ma- 
trices Pi, Ti, bi operating in the "band", Gor'kov-Nambu 
and spin spaces, respectively; po, fo, (Tq being the corre- 
sponding unit matrices. 

We assume that the functions ^1,2 (p) have the form 
Ci,2(p) = T ^(p) + (^Mj where the parameter bpi de- 
scribes the deviation from the ideal nesting depending 
on doping. Thus, the band 1 and 2 are the hole and 
electron bands, respectively. The energy ^(p) can be lin- 
earized near the Fermi surface and we write it in the form 
C(p) = I'F (|p I —Pf) with = pj^/m assuming for sim- 
plicity that the Fermi velocities vp ~ |v| in the bands 
are equal to each other. Actually, this assumption is not 
fulfilled in real pnictidcs but the results obtained under 
this assumption are valid, at least qualitatively, even in 
the case of unequal Fermi velocities. 

For convenience we introduce new operators, da and 
6a, related to the operators cJjq,, c„a 



where 




(2.3) 



Then, one can write the Hamiltonian, Eq. (j2.ip . as 
follows. 

The kinetic energy part 1-ikin reads as 



Hfc™ = ^{(<7o)a/^[e(p)(4a/3 + 6^.^/3) (2.4) 
6fi{a'lap -bl^bp) + h.c.^ . 



The term describing the superconducting pairing takes 
the form 



nsc = -l {K^2)ap{^aaial + Abbibl-h.c.)} , 

(2.5) 

where Aa = A* and A^ = A2. 

The term related to the SDW can be written as follows 



T~isDW = 



[i{o2)o.pWMQaibl + h.c^ . (2.6) 



ba — C2a , 



(2.2) 



At last, the Zeeman term can be written as 



'Hz = -Y |w^Zl(((73)a;3C0s6'- ((Ti)a/3Sin6')djja/3 -hTyz2((0'3)Q^C0s6'-H ((Ti)a/3Sin6')fejj6^ (2.7) 
^ W^Z |(((73)a/3COs6'- (0-i)a/3Sin6')djja^ -f (((T3)a^COs6'+ ((Ti)a/3Sin6')feJ,6,3| , 



r 



where we set Wzi — Wz2 — Wz — HeffH. 

Eqs. I|2.4H2.7I) are cumbersome and it is not convenient 



to use them directly. Fortunately, they can be rewritten 
in a more compact form introducing the operators Ana 
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and Bna which are matrices in Gor'kov-Nambu (index n) 
and spin (index a) space (see for exampleS,): 



La — ; 



iaa = al ; (2. 
B2a - bi . 



In order to take into account two bands or, in other 
words, different parts of the Fermi surface, we define op- 
erators Cmna with the index m related to different bands 
so that 



(2.9) 



The operators Cmna obey the commutation relations 

C'Lc.C'jn^ + ClnpClma = ^klSmnSap , (2.10) 
CkrnaClnf} + ClnpCkma = ^klSmn^a^ ■ (^-H) 

After rewriting the energy terms Eqs. (12. 4112. 7p in 

terms of the operators {C'm„Q, Cmna} = {(7^,(7} the 
Hamiltonian in Eq. (j2.1[) can be written in the following 
way 



(2.12) 



Hfcin = ^(p) ■ Po ■ T3 ■ aa ~ Sfi ■ p3 ■ Ts ■ ao , 

jj^^ _ I A'po • fi • (73 + iA"p3 ■ iT2 ■ (73 , S++ 
I A'p3 • fi • (73 + iA"po ■ iT2 ■ (73 , s+_ 

Hz = -~Wz ■ {po • To • (73 COS 6* - p3 • fs • (7i sin 6*} 



(2.14) 
(2.15) 
(2.16) 
(2.17) 



The order parameter A = A' + lA" is related to Ai^2 
as Ai = A* = A2 (s++-pairing) and Ai = A* = — A2 
(s+_-pairing). 

One can perform a rotation in the spin space and 
change the direction of the field Hea, and the magnetiza- 
tion m. The rotation around the axis j — x,y, z by the 
angle i? means a unitary transformation of the Hamilto- 
nian H, 



H RjHr] , ^ (2.18) 
where the unitary rotation matrix Rj has the form 



and 



R-j = cos(?9/2) + fj sin(i?/2) 



I • P3 • T3 • CTi , J = x 

, J - X i • P3 • fa • (72 , j = 2/ 
I • po • To • (73 , J 



(2.19) 



(2.20) 



If the magnetization in the SDW is oriented in the x- 
or ^/-direction, then Hsdw has the form 



where the operator H has the form 



H = Hfe„ + Bsc + HsDW + Hz 



with 



(2.13) 



iisDW{x,y) — iW^MO • /52 • T2 • (7i,2 



(2.21) 



Having specified the Hamiltonian of the model in the 
compact matrix notation we can introduce the Green's 
functions in terms of the operators & ,C in the usual 
way. For example, the retarded Green's function can 
be written as 



G«(p, p'; t, t') = [l/i) (C(p; t) ■ C\p'; t') + C+(p'; t') ■ ^(p; t)) e{t - t') , 
while the Keldysh function takes the form 

G^(p, p'; t, t') = (1/z) (C(p; t) ■ C^p'; t') - C^p'; t') ■ (7(p; t)) . 



(2.22) 



(2.23) 



One can define the matrix Green's function G with the 
block elements G^, G^' and G^ (see^iiiii^-^) 



G = 



G" G^ 
G^ 



(2.24) 
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Using these functions we can calculate various macro- 
scopic quantities. The Green's functions obey dynamic 
equations analogous to the Eilenberger equation and 
their generalizations to the nonequilibrium casci^ i'^^'^^" — 



III. EQUATIONS FOR QUASICLASSICAL 
GREEN'S FUNCTIONS 



Following the method of the quasiclassical Green 
functions^i^i^""— one should introduce the quasiclassical 
matrix Green's function g related to the Green's function 
G(r, r; t, t') as 



{i/ir) / di (f3 • G) 



(3.1) 



and derive dynamic equations for it. This section is de- 
voted to such a derivation that can be carried out in the 
standard way. 

The equations for the matrix G can be obtained from 
the equation of motion for the operators C(k; t), 



(3.2) 



and the definition of the Green's functions G^, G^ and 
G^, Eqs. (|2.22H2.23l) . 

Taking into account the commutation relations (j2.1ip , 
we obtain 



idtG - H • G = I5{t - t') 



(3.3) 



where the matrix H is defined in Eqs. (|2.13fET7|) . 

When deriving Eq. p.3|) . we used the property 



H 



H 



Analogously, one can obtain the conjugate equation 



idt'G-G-Yl = l5{t-t') 



(3.4) 



The next steps are standard for deriving the Eilen- 
berger equatioit^ or the equation for a more general 
Green's function GML^^. We multiply Eq. by 
f3 from the left and Eq. p.4p from the right and sub- 
tract these equations from each other. Then, the ob- 
tained equation is integrated over the energy ^(p) where 
P = l(P + P')l/2- Finally, we obtain the equation for 
the matrix Green's function 



(rs • dtg + dt'g ■ fg) + vVg + ^ [A^'^ , g] = , (3.5) 



where 



A 



U) 



'sc 



^SDW 



(3.6) 



with ("±" stands for s_| -pairing) 



■ P3 ■ To ■ o-Q , 



(3.7) 



hfp = -A' • /5o,3 • ■ 0-3 - iA" • p^^ • n • (T3 , (3.. 



^SDW 



i ■ Wmo • /52 • n • (Ti , j = x 
= { ~i- Wmo ■ h ■ n ■ ^2 , j = y 



^ i ■ Wmo • Pi • T2 • cts , 



(3.9) 



(3.10) 



where the matrix P — Pi + P3 with Pi = ps ■ tq ■ ai sin 9 
and P3 = —po-T3-a3 cos 9 describes the external magnetic 
field. 

Using the Matsubara representation one reduces 
Eq. ((3?5|) to the form 



[uJnT3,g] +vVg + i[A'£ rg] =0. 



(3.11) 



Repeating arguments used in the derivation for 
superconductor a^'^^i'*^' '— or two-band metals^ a normal- 
ization condition for the matrix g can easily be derived 



g{t,ti)-giti,t')^16{t-t') 



(3.12) 



with = S,nnSstSa0- 

Eqs. (|3. 5113. 12)) supplemented by proper boundary con- 
ditions allow one to find unambiguously the quasiclassical 
Green's function g. As soon as this matrix Green's func- 
tion is known, one can calculate macroscopic quantities 
of interest. For example the current density jcj in the 
system is equal to 



L = -^eiy J deTr{p3 • fg • (Tq • (v.g(e,e'))) , (3-13) 

where e = (e -I- ui — t ~ t' , v — pFin/i:'^ is the 

density-of-states at the Fermi energy and the angle brack- 
ets mean the averaging over the momentum directions: 



((■••))=/§(•■•) 



(3.14) 



The components of the static magnetic moment M in 
the a;-, y- and z-directions are given by 



^ 00 

M^^y ^ Mq ~ i{2nT)iypB-^ ^ Tr(/53 • fo ■ (Ti,2 ■ g(a;)) , 

uj=0 

(3.15) 

00 

= Mo - i{2TTT)i'pB-^ ^ Tr(/5o • f3 ■ 0-3 ■ g{uj)) , 

uj=0 

(3.16) 
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where Mq = y^^^jvH is the Pauh paramagnetic term. 
This term arises as a result of the integration over mo- 
menta far from the Fermi surface^ and cannot be cal- 
culated in the quasiclassical approximation. Writing 
Eqs. (|3.151 13.161) we replaced the integration over the 
energy e by the summation over the Matsubara frequen- 
cies. Below we are interested in the polarization of elec- 
tron spins by a static magnetic field H and this is why 
Eqs. (I3.15[ I3.16P are written for the static case. 

Using Eq. p.5|) one can obtain the expressions for the 
spin currents j^p. For example, the spin currents with 
(x, ?/)-spin projections are given by 

oo 

= -^{2^T)-^iBV Tr(p3 • h ■ ^la ' (v3(e, • 

aj=0 

(3.17) 

The order parameters are defined from the conven- 
tional self-consistency equations 

Ai,2 ^ / rfeTr((/5o ± ps) • 17-2 • 0-3 • (5(e,e'))^) , 
o Jo 

(3.18) 

Wmz = / rfeTr(pi • 1x2 ■ 0-3 • (g(e, e')) ) ■ 

o Jo 

(3.19) 

Eqs. (P3HXTU|) . together with Eqs. ((XTBHXTO)) and the 
expressions for the electric current, Eq. p. 131) . and mag- 
netization, Eqs. (|3.15H3.16)) . allow one to study various 
problems both in homogeneous and nonhomogeneous sys- 
tems. 

In the next section we calculate the magnetic moment 
induced in the system by an applied magnetic field and 



and uJn — TTT{2n + 1) is the Matsubara frequency. The 
matrix P is defined in Eq. p.lOp . 

For simplicity, we neglect the deviation from the per- 
fect nesting and set Sfj, = 0. We also consider the case 

of the S-i -pairing and the magnetization being oriented 

along the z-axis (j = z). The energy of the Zeeman split- 
ting is assumed to be small, Wz ^ {A, Wj\/o}, which 
allows us to consider the right-hand side of Eq. (|4.ip as 
a perturbation. 

In the zeroth order in Wz we neglect the R.H.S. of 
Eq. (j4.1l) and obtain the homogeneous equation 

[5o,M'^]=0. (4.3) 



find the Knight shift. 



IV. THE KNIGHT SHIFT 

In the experiments,'^^ the nuclear magnetic resonance 
(NMR) in superconducting (Sn) granules was studied. It 
was found that at low temperatures T the resonance line 
is shifted with respect to its position in the absence of the 
electron polarization (the so called Knight shift). Since 
at low T the free electrons in tin are bound in singlet 
Cooper pairs, they cannot contribute to the magnetic 
moment of the granules. Abrikosov and Gor'kov^ sug- 
gested an explanation for the observed Knight shift tak- 
ing into account the spin-orbit interaction. They showed 
that even at zero temperature this interaction gives rise 
to a non-zero polarization of electron spins in an external 
magnetic field H. 

In this section we study the Knight shift in a supercon- 
ductor with an SDW and show that even in the absence 
of the spin-orbit interaction the Knight shift is finite pro- 
vided the magnetic field H is not parallel to the direction 
of the magnetization in the SDW. 

In order to calculate the electron spin polarization in 
the field H we use Eq. p.5|) for the retarded (or advanced 
Green's functions) written in the Matsubara representa- 
tion, Eq. (|3.1ip . Since we consider the uniform case, the 
second term in Eq. p. lip may be omitted. Thus, we have 
to solve the equation 



[g,n^^]^[g,tWzP], (4.1) 

where 



= X 

= y (4.2) 

= z 

I 

In principle, any matrix function of fl^i^ satisfies 
Eq. (|4.3p . In order to find the function go unambigu- 
ously one should check whether the solution satisfies the 
normalization imposed by Eq. p.l2p or not. This leads 
us to the solution 



go = ^M'^ , (4.4) 

where £g — lu"^ + VK^fo + I^P- I* easy to see that the 
solution given by Eq. (|4.4p really satisfies the normaliza- 
tion condition 



I 

P2 • n • (72 , i 

-Pl • ■r2 • 0-3 , j 
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FIG. 1. (Color online.) Temperature dependence of the magnetic moment induced by a weak magnetic field for different ratios 
of magnetic to superconducting order parameters. 

a) M^T): ^ as T ^ 0; 

b) M^{T): stays finite as T -> 0. 



50 ■ 50 = i • (4.5) 
The correction Sg ~ g ~ ga obeys the equation 

[,55 , fl^'^] = £s[Sg , go] = -2£sgo ■ Sg = iWz [go , P] . 

(4.6) 



Here, we used the relation 



here those parts of 5g that contribute to the magnetic 
moments, i.e. Sgx and Sgz'- 



lAP 

Sgx = Sgi = -iWz-p:^h -tq-oi sm6' , (4.9) 



% = = -iWz ^ ' ^^° po • ■ ^3 COS0 . (4.10) 



Using Eqs. p.lSffOG]) we find for the spin polarization 
induced by the field H^x 



Sg ■ go + go ■ Sg ^0 , 

which follows from the normalization condition. 
Using Eq. (|4.5p we find from Eq. (|4.6p 



Sg^~'^{go-P-go~P) 



(4.7) 



(4.8) 



The matrix Sg can be written in an explicit form with 
the help of Eq. (14. 4p and the expression for P. We present 



lAI 



(4.11) 



LJ^O 



= Mo I 1 - 2TTT ^ +3^'^° I COS0 . (4.12) 



Ul^O 



In principle, the sums over u! in Eqs. (|4.11[HTT^ can be 
calculated at arbitrary temperatures but the final expres- 
sions are somewhat cumbersome. Therefore, we restrict 
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ourselves by the limit of low temperatures T — s- 0. In 
this limit one replaces the sums over uj by integrals, which 
leads to the following expressions 



As a result, we obtain 



sm f 



= Mo 1 



^1 



(4.13) 
) cos6' ^ 0. (4.14) 



Eq. (|4.14p shows that in the limit of the perfect nesting 
((5/i = 0) the Knight shift vanishes at T = provided 
the applied magnetic field H is parallel to the orientation 
of the magnetization of the SDW. If the direction of H 
deviates from the ^-direction, the Knight shift is finite 
and the induced magnetic moment approaches the 
spin magnetic moment of free electrons (Pauli paramag- 
netism) for Wmo 3> A. The obtained results do not 
depend on the relation between Ai and A2 and thus are 
valid for both s-j - and s++- pairing. 

In Fig. (HI) we plot the temperature dependence of 
and Mz for different ratios Wmq/ IS. assuming that in the 
considered temperature range the magnetization m of 
the SDW depends only weakly on T. 

At present, it is not easy to quantitatively compare 
our results obtained within a simplified model with avail- 
able experimental data concerning the NMR studies in 
pnictides^^i^. There are several reasons for the difficulty 
and the major one is that the Knight shift is actually 
not discussed in those papers. Furthermore, the influ- 
ence of the free electrons on the position of the NMR 
peak in the compound NaFeAs containing ^^Na atoms 
is weak because free electrons move in the FeAs planes. 
On the other hand, an internal local magnetic fleld of 
the SDW causes a stronger influence on the NMR peaks 
corresponding to the atoms ''^As. In addition, the ideal 
nesting is assumed in our model that leads to the fully 
gapped Fermi surface. This assumption is not fulfilled in 
materials studied in Refsi^i^. 



PROXIMITY EFFECT 



(5.1) 



where rix = Vx/v = cos a, a is the angle between the 
Fermi velocity v and the x-axis, and v = |v|. The 
matrix II is specified in Eq. (|4.2p . 

We assume the s_| -pairing and let the magnetization 

be directed along the z-axis. For simplicity, we assume 
as previously perfect nesting by putting 5/i — 0. 

Two different cases will be considered now: 

a) A'(x) = Ao, A"(a;) = 0, Wmo{x) = Wmo,s for 
a; < and A'(x) = A"(a;) = 0, Wmo{x) = Wmo,m 
for a; > 

b) A'(x) = Ao, A"ix) = 0, WMoix) = Wmo,s for 
a; < and A' (a;) = A"(x) = 0, Wmo{x) = for 
a; > 

The second case corresponds to an interface between a 
superconductor with SDW and a normal (nonmagnetic) 
metal. We denote this type of contacts as Ssdw/N. 

The first case corresponds to a system with SDW hav- 
ing the superconducting order parameter at a; < 0. This 
type of contact is denoted as Ssdw/^sdw- 

The contact between the two regions is assumed to be 
ideal and therefore all the functions should be continuous 
across the boundary a; = 0. 

As in the previous section, we represent the solution 
^(a;) in the form 



(5.2) 



Here the matrix qq is a constant in space and obeys the 
equation 



[50, M^^] =0 



(5.3) 



The proper solution of this equation is written again in 
the form 



We study the proximity effect considering a simple 
model: a contact between a superconductor with SDW 
(two order parameters: A and Wm), which we denote as 
SsDW, and a conductor (or insulator at Wm 7^ and 
low temperatures) with one order parameter (for example 
VFm) or with a simple normal metal iV (A = Wm = 0). 
Such a case may be realized in pnictides with a nonuni- 
form doping level. We will find the quasiclassical retarded 
and advanced Green's functions describing the equilib- 
rium properties such as the density-of-states or the order 
parameters, A and Wm- 

These Green's functions obey the generalized Eilen- 
berger equation. This equation is obtained from Eq. (13. 5p 
by taking its element (11) or (22) and performing the 
Fourier transformation in the Matsubara representation. 



90 = -^n 



(z) 



(5.4) 



with £g = UJ, 



2 I A2 -I- U^2 



The matrix 6g{x) is not supposed to be small. It can 
be split into an even, s, and odd, d, in parts 



dg ^ s + 71x0, . 



(5.5) 



Substituting Eq. (|5.5p into Eq. (|5.1I) and separating 
the even and odd in parts, we come to equations 



vdxs + [fl^^^ ,a] = , 
nlvdxci + [ll!f^ , s] = . 



(5.6) 
(5.7) 
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One can exclude the anisotropic part a by differenti- 
ating Eq. (I5.6P with respect to the coordinate x. Then, 
the equation for the isotropic part has the form 



2 9^ 



2£l 



Els 



= 0, 



(5.8) 



where laS = nxv/2£s is a characteristic length of pen- 
etration of perturbations caused by the proximity effect 
into the superconductor with SDW. At low temperatures 
this length is determined by the smallest of the lengths 

{^A =i w/Ao,Cm =i v/Wmo.m}- 

In the region a; > the characteristic length is 
laM = nxv/2£M (the case a)) or I^n = n^vj^E^ (the 
case b)), where E\i = -I- VF^fCM ^.nd = w^^j. 

We look for a solution of Eq. (|5.8p in the form 



I A'\'^)h + B{x)p-i • fa • (73 C{x)pi • fa • (73 , a; < 

yA(x)T-i -f B{x)p3 • f 2 • (73 -I- C(x)pi • f2 • (73 , X > 

(5.9) 

Substituting Eq. (|5.9p into Eq. (|5.8p one can rather 
easily find the functions A, B, C, etc. 



A{x) 


= Ao exp(x//Q5) , 


A{x) — exp(- 


-x/lau) , 

(5.10) 


B{x) 


= Bo exp(a;/Zas) , 


i?(x) — Bo exp(- 


-x/lau) , 

(5.11) 


C{x) 


= Coexp(a;/?as) , 


C(x) = Co exp(- 


-x/laM) , 

(5.12) 




AoWn -1- BoAo - 


1- CoWmo,s = , 


(5.13) 






1- CoWmo,s = . 


(5.14) 



We see that there are four independent arbitrary con- 
stants: Bq, Co, Bq and Co characterizing the solution 
s. They should be determined from the matching condi- 
tions that require the continuity of the matrices go+s and 
d. The matching conditions are reduced to six equations 
four of which are independent and the other two follow 
from these four equations. Solving these equations, we 
find 



Ao 



Ba = - 
Co 



i;mM - ms 



Aq^ — T— — , (5.15) 



1 + T 



1 + ^ 



So = 



-A, 







l + <r 
<,ms - niM 

1 + T 



(5.16) 
(5.17) 



where ^ = (wst^m + msuiM) ^, t^s = ^/Es, 
zum — lu/ Em are the normalized Matsubara frequencies 
in the SC and M regions, and ms.M — Wmqs.m /Es,m 
are the magnetic order parameters in these regions. In 
the case of a contact of a superconductor with SDW and 
of a normal metal {Ssdw/N contact), the energy Em 
should be replaced by Eu and the quantity Wmo set to 
be zero. 

The amplitudes Bq, Cq determine the penetration of 
the superconducting and magnetic correlations into the 
region with SDW or into the normal metal N due to 
the proximity effect. The amplitudes Bq, Cq describe 
the inverse proximity effect or, in other words, a sup- 
pression of Ao and Wmo in the superconductor near the 
SsDW /NsDW interface (or in the Ssdw/N interface) 
due to the inverse proximity effect. 

Note that, strictly speaking, in the case of the 
SsDW /NsDW system we have to calculate the mag- 
netic order parameter Mmo self-consistently using the 
amplitude Cq. This makes the problem more difficult. 
However, the obtained results remain valid provided the 
quantity Wmo is the same at a; > and a; < (i.e. 
Wmo,S = Wmo.m = Wmo)- In this case we obtain 



An = - 



Bn^- 



Co 



EsEm{Es + Em) 

^oEm 
Es{Es + Em) ' 
moA^ 
EsEm{Es + Em) 



io = 0, (5.18) 
Ao 



So- 



(5.19) 



Es + Em 
Co = . (5.20) 



The results obtained mean that the corrections to the 
DOS and to the magnetic order parameter determined by 
Aq and Co are absent in this case. The superconducting 
pair function B{x) penetrates the region with SDW over 
the length laM with the amplitude Bq. As it should be, 
this pair function penetrates the N region over the length 

laN = lau = \nx\v/uJn- 



VI. JOSEPHSON EFFECT 

In this section, we calculate the dc Josephson cur- 
rent in an Ssdw/N/ Ssdw system using a simple model. 
We assume again an ideal nesting {S/j, = 0) and take 
into account the impurity scattering in the self-consistent 
Born approximation. Then, the Eilenberger equation, 
Eq. (j5.ip . for the matrix g acquires the form 



(6.1) 



where the angle brackets stand for the average over the 
directions of the momentum and the matrix 77 in the 
right (left) superconductors is equal to 
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n 



, P0,3 ■ T2 • 0-3 COS 



(!) 



± hfi ■ n ■ CTs sin 



in(|)) +M^Mo( 



MQ Pi • T2 • 0-3 COS 



(f) 



=F P2 ■ n ■ 0-2 sm 



■"(I): 



(6.2) 



N 



(6.4) 



where 




FIG. 2. (Color online.) Josephson junction under consider- 
ation. A normal metal film is placed between two supercon- 
ducting leads with SDW. The orientation of the magnetiza- 
tion vectors of the SDW in the leads lying in the y-z-plane is 
shown in red. 



Eq. (|6.2p corresponds to the case when the supercon- 
ducting phase in the right (left) Ssdw equals ±1^9/2 and 
the angle between the vector of the magnetization of the 
SDW and z-axis is equal to ±a/2 (cf. Fig. ©). Consid- 
ering the impurity scattering, we neglect interband scat- 
tering and regard, for simplicity, the impurity scattering 
time r equal for each band. In the middle of the iV layer 
A = and Wmo — 0. We assume that the scattering 
time T in this layer is rather short: hr"^ 3> T (diffusive 
limit). 

As a boundary condition, we adopt the one obtained 
for a simplified model suggested in Refi^ 



dx 



±7s [g , ~gs\ 



x=±L 



(6.3) 



x=±L 



where the parameter 7^ = (2i?B0'Ar)^^ depends on the 
interface resistance Rb per unit area and the conductiv- 
ity um of the normal metal. This parameter is assumed 
to be a scalar, which means that we neglect transitions 
between different bands at the interfaces. In a more gen- 
eral case 7b is a matrix^Sr.^. 

We assume that the proximity effect is weak, which 
corresponds to small values of the parameter jb£.n, 
where S^n — \/ D 12-kT. In this case the Green's func- 
tions in the Ss dw are only weakly perturbed by the con- 
tact with the N layer (because of a large interface resis- 
tance) and the first term in Eq. (|6.ip can be neglected. 
Then, the solution of Eq. (|6.1|) on the boundaries between 
the normal metal and superconductors can be written as 



9+ 



UJnh + Apo,3 • T2 • (73 COS (^^^ 

+ WmoPi ■f2- as COS j 



(6.5) 



and 



Ap3,o • n • (73 sin (^^"j (6.6) 
- WMap2 • n • (72 sin (^^^ 



with £l — ojI 



To find the Josephson current jj, we have to solve 
Eq. (pH in the N layer, where A = and Wmo = 0. 
As usually, we represent the Green's function g as 
a sum of the symmetric and antisymmetric functions: 
g =^ s + UxO, (see Eq. (j5.5p ). For s and a we obtain 
from Eq. (|6.ip the following equations 



vdxS + UJn[TS ,a] 

v{nldxd) + ujn [t3 , (s)] = , 



1 

'2^ 



{sd — as) , (6.7) 
(6.8) 



where the angle brackets denote the angle averaging. In 
the diffusive limit considered here {tT/H ^ 1) the sec- 
ond term in Eq. (16. 7p in the left hand side is small. The 
right-hand side can be transformed using the equation 



sd + as — , 



(6.9) 



which follows immediately from the normalization condi- 
tion Eq. (|3.12l) . Then, we obtain for d 



d = —Is dr 



(6.10) 



We used another part of the normalization condition (for 
symmetric function) 



-2 I 2-2 1 
s + n„a ~ 1 



(6.11) 



in which the second term can be neglected (as follows 
from Eq. (I6.1ip \d\ <^ 1). In the considered limit of 
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a weak proximity effect (|s| ^ 1) the matrix s can be 
represented in the form 



s = fa sgn(a;) + 5s , (6-12) 
where the matrix Ss « (Ss) obeys the equation 



d^Ss 



+ K^,5s = 



(6.13) 



with = lujnjDi D = vl/3 being the diffusion coef- 
ficient which is assumed to be the same in each band. 
Eq. (|6.13p follows from Eq. (|6.8I) because, as we will see, 
the matrix Ss anticommutes with the matrix fa. 

Eq. (|6.13l) has to be solved with the boundary condi- 
tions which follow from Eq. (|6.3p and the representation 
Eq. (lOa 



dSs 
dx 



i±L) = ±7£ 



gs{±L) - fg 



(6.14) 



The solution of Eq. ()6.14|) can easily be found in the form 



Ss{x) 
where 



A- 



cosh(K^2;) 
sinh(KtjL) 



B ■ 



smhlK^x) 
cosh(Kcji) ' 



3.15) 



jj = jc sin ip , (6.20) 
where the critical Josephson current jc is given by 



jc = a^%L{2TTT) 



(6.21) 



where 6^ — k^L. This formula differs from the expres- 
sion for jc in an S/N/S junction only by the term W^^q 
in the energy Ss- The presence of this term leads to a 
suppression of the current jc provided Wmo is not small 
compared to the superconducting energy gap A. Note 
that, in our simple model, the Josephson critical current 
does not depend on the mutual orientation of the vectors 
of the magnetization in the left and right superconduc- 
tors SsDW- 

It is of interest to calculate the spin current in the 
Josephson junction. The spin current is given by 
Eq. p.l7p . Using Eq. (|6.10p one can write the spin cur- 
rent between the superconductors as 



4^1 = -^^^^^^ETr{p3-r3-^,.<5^9.5s} , (6.22) 



where the upper index means the spin orientation, and 
the lower stands for the direction of the current. Substi- 
tuting Eqs. (pl^MT7|) . we obtain from Eq. 



and 



B 



IB 

£s 



IB 



A/5o,3 ■ f2 • (73 cos 



(f) 



WawPi ■ T2 ■ 0-3 COS ( - 



A/53,0 • f 1 • (73 sin 



(f) 



WmqP2 • T2 • (72 sin ( - 



a 



(6.16) 



3.17) 



The Josephson current jj can be calculated using 
Eq. p. 131) and the expression for a (|6.10p . We are inter- 
ested in the part of a which contributes to the current. 
This part can be written in the main approximation as 



J sp 



, sinh(26'^) 



sma . 



(6.23) 

Eq. (|6.23p resembles the Josephson expression for the 
supercurrent. Both formulas contain sine of an angle. In 
the conventional Josephson formula this angle is equal 
to the difference of the phases of the superconductors, 
whereas the angle a in the expression for the spin current 
in Eq. (|6.23p determines the mutual orientation of the 
SDW in the right and left electrodes. Eq. (|6.22p is valid 
also above the superconducting transition temperature 
Tc, when A = 0. This current is dissipationless like the 
Josephson current. The nature of a similar dissipationless 
spin current in systems which differ from ours has been 
discussed in Refi^ 



d = —Ids dx5s . 



(6.18) 



Proceeding in this way we reduce the expression for 
the Josephson current jj to the form 



J J = ^ Tr{p3 ■ f3 • ,70 • 5s djs] . (6.19) 



Substituting Eqs. ((OBMT7|) into Eq. (I5T^ . we find 



VII. DISCUSSION 

We have derived equations for the quasiclassical 
Green's functions for a two-band superconductor with 
an SDW. It was assumed that the Fermi velocities in 
the electron and hole bands are equal. We neglected the 
anisotropy of the Fermi surfaces. 

Using these equations and assuming the ideal nesting, 
we considered three problems: the Knight shift, the prox- 
imity effect and the dc Josephson effect. It was shown 
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that, provided the direction of the applied magnetic field 
coincides with the direction of the magnetization m in 
the SDW, the Knight shift vanishes at zero temperature 
and in the absence of the spin-orbit interaction. If the 
magnetic field is not collinear with the m vector, the 
Knight shift is finite, which correlates with the results 
of a recent paper, where it was shown that the DOS 
of a superconductor with the SDW also depends on the 
orientation of the external magnetic field. 

We have demonstrated that near the interface between 
the superconductor with the SDW and a normal metal 
the components of g describing both magnetic and super- 
conducting correlations penetrate into the normal metal 
over the length of the order min{hv /Wmo, ftw/A}. 

Using the simplest model of the Ssdw /N / Ssdw 
Josephson junction, we calculated the critical Josephson 
current and showed that, in this model, it does not de- 
pend on the mutual orientation of the magnetization in 
the superconductors with the SDW. 

However, the dissipationless spin current (jsp), which 
arises in the junction, depends on the misorientation an- 
gle between the magnetisations of the SDW, a, in the 
same way {jsp ~ sin a) as the Josephson current jj de- 
pends on the phase difference of the superconducting or- 
der parameter. 

Although the equations for the quasiclassical Green's 
functions have been derived using the simplest model of a 
two-band superconductor with an SDW, we believe that 
the results obtained on the basis of this model remain 



valid, at least qualitatively, for more complicated models 
describing realistic materials. 

The derived equations can be easily generalized to the 
case of impurity scattering and can be applied to different 
problems, equilibrium and nonequilibrium (ac Josephson 
effects in Ssdw /N j Ssdw junctions, vortices, etc.). 

Note added in proof. Although Eqs. (|3.51 13.12p for the 
quasiclassical Green's functions for two-band supercon- 
ductors with the SDW are valid for arbitrary deviation 
from the ideal nesting ^ 0), we assumed the ideal 
nesting (/i — 0) when applying these equations to the 
study of particular effects. This assumption is justified, 
strictly speaking, in a hypothetical case of equal critical 
transition temperatures {Tcs — Tcsdw) or in the case 
of a metastable states (when a first-order transition takes 
place). In real materials 7^ 0, and therefore one has to 
take into account the dependence of the order parameters 
(A and Wm) and other quantities on /i. For example, one 
can show that the Josephson critical current jc contains 
an additional term 5jc which is negative and depends on 
the angle a: Sjc ^ — /i^ cos^ a. 
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